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Functionalization of nanoparticles or colloids is increasingly being 
used to develop customizable "atoms". Functionalization by at- 
taching single strands of DNA allows for direct control of the bind- 
ing between nanoparticles, since hybridization of double strands will 
only occur if the base sequences of single strands are complemen- 
tary. Nanoparticles and colloids functionalized by four single strands 
of DNA can be thought of as designed analogs to tetrahedral net- 
work forming atoms and molecules, with a difference that the at- 
tached DNA strands allow for control of the length scale of bonding 
relative to the core size. We explore the behavior of an experimen- 
tally realized model for a nanoparticles functionalized by four single 
strands of DNA (a tetramer), and show that this single-component 
model exhibits a rich phase diagram with at least three critical points 
and four thermodynamically distinct amorphous phases. We demon- 
strate that the additional critical points are part of the Ising univer- 
sality class, like the ordinary liquid-gas critical point. The dense 
phases consist of a hierarchy of interpenetrating networks, reminis- 
cent of a woven cloth. Thus bonding specificity of DNA provides 
an effective route to generate new nano- networked materials with 
polyamorphic behavior. The concept of network interpenetration 
helps to explain the generation of multiple liquid phases in single- 
component systems, suggested to occur in some atomic and molec- 
ular network forming fluids, including water and silica. 
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Technological developments are creating a wide variety of 
new building blocks which play the role of "functionalized 
atoms" for designed materials [1,2]. Borrowing ideas from bi- 
ological specificity [3-6] , it is now possible to control bonding 
between engineered building blocks in a selective way, mov- 
ing in the direction of an effective synthetic bottom-up strat- 
egy for self-assembly. The complementary binding of base 
pairs, combined with the ability to directly control the base 
sequence, makes DNA an ideal candidate for the development 
of network-based, nanostructured materials [5-8]. A variety 
of experiments [7-15] have demonstrated the possibility to re- 
alize nanostructured, DNA-based materials, including ordered 
crystal structures [9, 10]. 

By grafting short (i.e. smaller than the persistence length) 
single strands (ss) of DNA to core particles, the core can 
act as a "node" of a complex network formed when single 
strands combine into double-stranded (ds) DNA. Longer link- 
ing "arms" may be achieved by using a double-stranded spacer 
near the core nanoparticle [9, 10]. With an intelligent choice 
of the core, it is possible to control the nearest neighbor coor- 
dination of the functionalized nanoparticle, providing a direct 
route to study materials where the bonding coordination is 
much less than that of spherically symmetric molecules [16]. 
In such limited coordination systems, the general features of 
network formation are expected to be observed. 



By functionalizing the central core with four DNA strands, 
we expect to generate an engineered analog of tetrahedral net- 
work forming atoms or molecules, in which the bonding con- 
tribution to the interparticle interaction energy can be tuned 
by modulating the length L of the DNA strand, as well as 
the chemical properties of the solvent. Tetrahedrally coordi- 
nated fluids, such as water [17-21], silicon [22], and silica [23] 
appear to have unusual phase diagrams that include two crit- 
ical points: (i) the usual gas- liquid critical point, and (ii) a 
lower temperature critical point terminating the coexistence 
of a low-density and high-density liquid. This phenomenon is 
also referred to as "polyamorphism" [24]. Our model DNA 
"molecules" allow us to explore if such unusual polyamorphic 
behavior is shared by designed systems, and also provide in- 
sights that might help understand complex phase diagrams in 
traditional materials. Our findings indicate that interpenetra- 
tion of repeating, amorphous networks can give rise to multi- 
ple fluid phases, providing a disordered analog of the natural 




Fig. 1. Visualization of the model before and after coarse graining. A single den- 
drimer is formed by joining four single strands to a central hub (the four red spheres). 
The small green spheres are the attractive "sticky" spots that carry information about 
base type. The top panels show (A) a bonded pair of dendrimers and (B) a nearby pair 
that is unbonded. In the coarse-grained model (bottom panels), the 4 central sites are 
replaced by a single core (indicated by a blue sphere), and the arms (indicated by a 
blue cylinder), each representing a single strand of DNA, the arms are parameterized 
by a distance and angle dependent bonding potential. The schematics show both 
(C) bound and (D) unbound cases of two nearby dendrimers. The potential between 
dendrimers depends on the core-core separation r and the orientations Oi and O2 of 
the binding arms relative to r. 
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Fig. 2. Dendrimer-dendrimer pair potential for the model studied. (A) The po- 
tential (5V{r, 61,62 = 0) = (3{VNB{r) + (r, 6>i , 6*2 = 0)) along the 
ideal bonding direction, where (3 = l/ksT. For simplicity of the figure, we 
take Pb = 1, so that the full effect of the bonding potential is apparent. (B) 
The angular variation of the potential at the ideal bonding distance. The attrac- 
tion drops rapidly as the 61 departs from zero. Note that since labels are arbitrary, 
V{r = d,6i = 0,62) vs. 62 is identical to V{r = d, 61,62 = 0) vs 61. 

example of interpenetration which occurs in ices VI, VII, and 
VIII [25] , and which helps to understand how interpenetration 
in liquid water is feasible at high pressures. 

Our results are based on a series of Monte Carlo simula- 
tions of a simple, coarse-grained molecular model for tetrahe- 
dral DNA dendrimers - that is, a core nanoparticle binding 
four ss-DNA with tetrahedral symmetry. Such tetramers have 
also been experimentally synthesized [11]. The model retains 
two important basic features: (i) base-pair selectivity and (ii) 
the limitation that each ss-DNA can bind to only one other ss- 
DNA to form a ds-DNA. The model omits chemical details and 
explicit solvent effects which would make calculations of the 
phase behavior infeasible. Each molecular unit of the model is 
composed of a tetrahedral hub tethering four identical DNA- 
like strands composed of L connected monomers. The detailed 
model (introduced in refs. [26, 27]) is shown for the case L = 8 
in Fig. 1 A and B. 

Here we focus mainly on a simplified model [28] that is 
directly derived from the more complex model. We represent 
each dendrimer by a central core with rigid arms - represent- 
ing the ss-DNA - that point to the vertices of tetrahedron 
(Fig. 1 C and D). The coarse-grained dendrimer-dendrimer 
interaction potential is a function only of (i) the separation r 
of the cores of the nanoparticles, (ii) the relative orientations 
Oi and O2 of the binding strands, and (iii) the temperature T 
dependence of the probability Pb{T) that single strands bond 
to form double strands [28]. Fig. 2 shows both the separa- 
tion dependence of the potential along the bonding direction 
(^1 =^2 = 0) and the angular dependence at the bonding dis- 
tance. The details of the potential and the parameterization 
from the more detailed model are described in the Methods 
section. In contrast with atomic or molecular network form- 
ing systems, the core-core excluded volume interactions of our 
model are relevant only at separations much smaller than the 




Fig. 3. The distribution of the density P(p) for L = 4 at the locations (Tc, 
/ic) of the three distinct critical points terminating first-order fluid-fluid transitions. 
In the inset, we identify the universality class by calculating the distribution of the 
order parameter P(M), which coincides with the known distribution for the Ising 
model (solid black line). Note that M is shifted by its mean value and scaled so that 
it has unit variance. 

bonding distance, allowing for significant interpenetration ef- 
fects. The model was originally parameterized for DNA se- 
quences of length L = 8 [28] ; here we generalize the model to 
shorter (L = 4) and longer (L = 16) strands. The sequence of 
the DNA bases in each arm is chosen so that the complemen- 
tary sequence is identical, but in reverse order. As a result 
head-to-tail bonding between the arms of different dendrimers 
occurs, and when many dendrimers bind, they form an open 
tetrahedral network. 

Results 

To evaluate the phase behavior, we first determine the loca- 
tion of any critical points in the phase diagram. In the grand 
canonical ensemble, we can accurately calculate the position 
of a critical point by locating the state points (defined by 
chemical potential /i and T) where the density fluctuations 
exhibit a specific bimodal distribution. We define the density 
p = N/VdP , where d ^ L is the most probable separation of 
nearest neighbor dendrimer cores when they are bonded, N 
is the number of dendrimers and V is the volume. Including 
d in the definition of p simplifies comparison of systems with 
different L. Surprisingly, as shown in Fig. 3, we discover three 
different critical points for all investigated L, namely: the ex- 
pected critical point terminating the transition between the 
unassociated dendrimers and a networked fluid state (the ana- 
log of the gas-liquid critical point in fluids) as well as two ad- 
ditional liquid-liquid critical points. To demonstrate that all 
three critical points belong to the Ising universality class, we 
show that the calculated order parameter distribution P(M), 
where M — p — su, coincides with the known distribution of 
the magnetization in the Ising model at the critical point [29] . 
Here, u is the potential energy density and s is referred to as 
the field-mixing parameter. The possibility of multiple critical 
points in networked liquids is well documented [17-23, 30, 31], 
but it has never before been possible to quantitatively observe 
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Fig. 4. Phase diagrams and schematics of fluid structure. We plot the phase 
diagrams for arms of length (A) L = 4, (B) L = 8 and (C) L = 16 bases. 
The red circles indicate the critical point locations. The black symbols indicate the 
phase boundary determined by a histogram reweighting procedure [39], and the blue 
symbols are determined from grand canonical ensemble simulations. We illustrate the 
structure of the dense fluids in panels (D-F) which correspond the the low, medium, 
and high density networked states, as indicated by the arrows. Green cylinders are 
drawn between the centers of bonded dendrimers. The arms of the dendrimers are 
indicated only by a small white sphere to provide the molecular orientation. 

the critical behavior and identify the universaUty class of any 
additional transitions. Additionally, the data in fig. 3 unam- 
biguously show that there may be a hierarchy of many such 
critical points. 

To complete the picture of the phase behavior, we iden- 
tify the coexisting phases and evaluate the phase boundaries 
(Fig. 4 A-C). For each L, three adjacent phase coexistence 
regions emerge, with roughly same width in p, suggesting the 
presence of well-defined preferential p values for the stable 
dense phases. The p dependence for all L can be better un- 
derstood by scaling p by the corresponding ideal diamond lat- 
tice density pd = 3a/3/8 ~ 0.650 (in units where the nearest 
neighbor lattice sites are separated by unit length) [32]. In 
this representation, the densities of the dense stable phases 
are approximatively integer multiples of pd, suggesting that 
these phases consist of interpenetrating four-coordinated net- 
works. For sufficiently large L, where the core size is very 
small compared to the bonding distance, we expect the for- 
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mat ion of additional interpenetrating networks, and of their 
associated critical points, beyond those we have already ob- 
served. So far our studies have not identified such a scenario, 
due to the excessive computations required to equilibrate at 
large densities and small temperatures. Unlike density, scaling 
of the temperature to a universal form is more complicated. 
Indeed, the differences in the L dependence of the bonding 
entropy causes the relative value of the critical temperature 
Tc to differ. A graphical representation of the dense stable 
phases is shown in Fig. 4 D-F. 

To quantify how interpenetrating networks are manifested 
in typical measured structural properties, and to provide ev- 
idence for the absence of crystal order, we evaluate the ra- 
dial distribution g{r) (where r is the separation between the 
centers of dendrimers) and its Fourier transform, the struc- 
ture factor S{q) (Fig. 5 A and B). The tetrahedral order of 
the lowest density networked phase is evidenced by the fact 
that the ratio of the position of the first and second neighbor 
peaks of g(r) is approximately 4/ a/6 ~ 1.63, expected for ideal 
tetrahedral coordination. The second and third interpenetrat- 
ing liquids also have the same ratio between first and second 
neighbor peak locations, but there is increased amplitude of 
g{r) for r/d < 1. Such an increase is consistent with the 
presence of dendrimers filling the holes of the complementary 
networks. The presence of neighbors closer than the bonding 
distance is made possible by the absence of significant hard- 
core interactions. Similar tetrahedral signatures are indicated 
by S{q). For the single- networked fluid, the correlation among 
the holes of the network gives rise to the "prepeak" in S{q) 
at q/qo ~ a/6/4, where is the location of the main peak 
of S{q). As these holes are filled by subsequent interpene- 
trating networks, the prepeak weakens, making it difficult to 
recognize from S{q) the underlying tetrahedral structure. 

To demonstrate that even the network of the densest phase 
is locally tetrahedral and interpenetrating, we look at the con- 
tribution to the g{r) restricted to dendrimers that are sepa- 
rated by a specific number of bonds, refereed to as chemical 
distance D (Fig. 5 C). More specifically, g{r) for D = 1 
only shows correlations with the nearest connected neighbors, 
while D = 2 includes both first and second connected neigh- 
bors, etc. In this way, we consider only neighboring particles 
that both belong to the same bond network as the selected 
central particle and are separated by no more than D bonds. 
The ordinary g{r) is recovered in the limit D ^ oo. Fig. 5 C 
shows that for D < 3, the system has nearly perfect tetrahe- 
dral order and there is almost no interpenetration. For D > 3, 
connections between networks destroy the long range order 
and dendrimers can loop back to very small physical distance 
r, giving rise to the increase in the amplitude at r/d < 1. We 
visualize the interpenetrating structure of the three networks 
in the densest liquid by rendering dendrimers and the bonds 
between dendrimers of distinct networks with distinct colors 
(Fig. 5 D-F). 

Lastly, we consider the possible concern that additional 
phase transitions may be only an artifact of the simplified 
model. Since comprehensive phase diagram calculations for 
the original model of fig. 1 (A,B) are infeasible, we instead 
check for phase separation in both the original and effective 
potentials by simulating at fixed total p in the phase coex- 
istence region. Fig. 6 shows S{q) after a long equilibration 
to allow phase separation to complete for both the original 
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Fig. 5. Demonstration of the existence of interpenetrating networks and quantification of their influence on traditional structural measures. (A) The radial distribution g{r) 
as a function of scaled distance r jd for three distinct liquids phases. The increased amplitude for r j d <i 1 is a consequence of the interpenetration. (B) The structure factor 
S'(g) for the same three liquids. (C) g{r) restricted to bonded neighbors up to a specific chemical distance D, providing evidence for the tetrahedrality of the interpenetrating 
networks. (D-F) Visualization of the state point used for calculations of c to show the three interpenetrating networks. Specifically, for an arbitrary selected dendrimer and for 
the two closest unbonded neighbors of that dendrimer, we draw their bonded neighbors up to D = 3. 



and effective potentials in the liquid-gas coexistence region, as 
well as the first liquid-liquid coexistence region. The growth 
of S((i) for small q is an indication of phase separation. The 
growth is much more pronounced for the liquid-gas coexis- 
tence because the ratio of densities of the coexisting phases 
is much larger than for the liquid-liquid coexistence. The be- 
havior of S{(i) for the real and effective potentials is nearly 
quantitatively identical, indicating that the same phase sepa- 
ration process occurs in the original potential. 



Summary 

The emergence of multiple critical points in this system re- 
sults from the possibility of creating multiple interpenetrat- 
ing bonded structures, explaining why such phenomenon is 
not observed in simple atomic and molecular systems for 
which hard-core interactions prevent significant interpenetra- 
tion. However, interpenetration is in principle possible in 
some networked liquids where directional bonding favors lo- 
cally "empty" bonded states, like water. In this respect, our 
work supports the possibility of multiple liquid-liquid criti- 
cal points in water and other networked liquids as a result of 
interpenetration of distorted bonded networks. Some care 
must be taken in applying interpenetration concepts to wa- 
ter and other molecular systems. Specifically, the hydrogen 
bond distance in water is not substantially larger than the 



core (oxygen-oxygen) repulsion. Despite this repulsion, wa- 
ter is able to arrange, at least in crystal forms VI, VII, and 
VIII, as two interpenetrating fully bonded lattices [25]. Un- 
like our DNA functionalized particles, the absence of bond- 
selectivity in water molecules means that molecules can share 
their bonding sites to form simultaneously pairs of bonds [33] . 
These shared bonds are absent for DNA functionalized parti- 
cles, where the lock-and-key mechanism imposes a one-to-one 
correspondence between bonding sites and number of bonded 
neighbors. Such shared bonds in water give rise to significant 
distortions of the network, which is clearly indicated from 
structural measurements [34]. As a result, in the case of wa- 
ter, the density of the high density interpenetrating phase 
should be significantly smaller than twice that of a single 
networked phase. Indeed, the difference between the low- 
and high-density states is experimentally known to be around 
25% [35]. Another approach that has helped to understand 
multiple phases in these dense systems is the use of spher- 
ically symmetric models with two length scales [19,36-38]. 
Here we have shown that a single bonding length scale that 
allows for interpenetration is sufficient to understand multi- 
ple amorphous phases, providing a complementary approach. 
Moreover, our results demonstrate that novel nanostructured 
materials can serve as model systems to help understand the 
behavior of conventional materials. 
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Fig. 6. Comparison of the structure for the original and effective potentials in the 
phase coexistence region for L = A. We calculate S{q) only after a long equili- 
bration to allow the phase separation process to near completion. The solid line is 
the effective potential, and the dotted line is the original potential of refs. [26,27]. 
The system density is fixed so that the system phase separates into coexisting phases. 
For the liquid-gas coexistence, we study T = 0.083 and p/ pd = 0.58; for the 
liquid-liquid coexistence, we study T = 0.079 and p/pd = 1-42, which can be 
readily identified in fig. 4 A. 

In summary, the chemical synthesis of new functionahzed 
nanoparticles offers the possibihty to control the bonding ori- 
entation and to generate bonding distances which are signif- 
icantly larger than the core size, a key element to favor the 
formation of multiple interpenetrating networks at high den- 
sity. DNA functionahzed cores provide a very effective way to 
realize such nanoparticle building blocks and to create struc- 
tured materials with properties that do not ordinarily appear 
in nature. Thus, our work suggests that it is possible to ex- 
pand the tool-box of building blocks so as to open the possi- 
bility of a hierarchy of networked phases, which could poten- 
tially be amorphous or ordered. In the present case, the inter- 
penetrating networks are identical. However, by appropriate 
mixing of particles functionahzed with distinct sequences of 
DNA, it should be possible to generate interpenetration in 
which each network has distinct properties. Such a material 
would provide three-dimensional nanoscopic canvass, woven 
by differently "colored" wires, the properties of which can be 
controlled directly by base sequence selection. 

Methods 

We first briefly describe the model and procedure for obtain- 
ing the parameters of the effective potential we use. The po- 
tential we study is a directly derived from a more detailed 
model [26, 27] for four-armed DNA dendrimers. Obtaining 
potential parameters for our coarse-grained representation re- 
quires an initial study of high and low T states in the more 
complex model. The more complex model includes explicit 
bases, and we choose a palindromic sequence of bases to en- 
sure head-to-tail bonding that gives rise to large networks. 
The bases comprise of four types of particles labeled A, C, G, 



and T. Bases of type A bond only with type T, and type C 
bonds only with type G. Starting from the tetrahedral core, 
the specific ordering of the bases used to determine the pa- 
rameters of the effective potential are: (i) L — 4, A-C-G-T; 
(ii) L = 8, A-C-G-T-A-C-G-T; (iii) L = 16, A-A-G-C-C-A-G- 
T-A-C-T-G-G-C-T-T. 

The coarse graining procedure has been explained in detail 
in Ref. [28]. We express the effective pair potential as sum 
of the potential for the bonding state Vb (weighted by the 
probability pb that two strands are bonded at temperature 
T) and the potential describing the non-bonding interaction 
Vnb, i-e. 

V{r, Oi,02,T) = Vnb (r) + pb (T) (r, , ^2) • [ 1 ] 

The potential Vnb is assumed to be spherically symmetric 
and it is very soft since it models primarily the repulsion of 
the small core. The potential Vb retains orientation depen- 
dence, since bonding only occurs when strands are properly 
aligned. To mimic bond selectivity and to reproduce the lock 
and key binding characteristic of biological binary associations 
like DNA, each arm may form no more than one bond. There- 
fore, if the arms labeled by Oi or O2 are already engaged in a 
bond with some other dendrimer in the system, Pb{T) = 0. 

To numerically evaluate the functional form of Vnb, Vb 
and Pb{T) we simulate a system composed of two dendrimers 
in a box for several T in the region where ss-DNA pairs. The 
numerical results for the probability Pb{T) can be modeled as 
a two-state expression 

where AU and AS measure the change in energy and entropy 
associated to the formation of a double strand and kB is the 
Boltzmann constant (in our reduced units, kB = !)• The an- 
gular dependence can be model as a Gaussian function. The 
ability of the effective potential to reproduce structural and 
thermodynamic properties of the explicit model has been ver- 
ified in ref. [28] for the case L = 8. We have confirmed that 
similar agreement in structural quantities is observed for the 
case L — 4 and L = 16, even in the region where multiple 
networks are observed. 

To evaluate the phase diagram, we carry out a series 
of Monte Carlo simulations in the grand canonical ensemble 
(fixed /i, volume V, and T) and study the distribution of the 
density fluctuations. To properly locate the critical points, 
we use histogram re- weighting techniques to ensure the order 
parameter distribution P[M) matches that of the Ising uni- 
versality class. The phase boundaries near the critical point 
are determined by a series of multi-canonical simulations, as 
described by ref. [39]. 
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